Anisotropic thermoelectric effect in helimagnetic tunnel junctions 
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Abstract 

Thermoelectric transport across normal-metal/helical- multiferroic/ferromagnetic heterojunc- 
tions is theoretically investigated. We find a anisotropic charge and spin thermopower with a Ci v 
symmetry. The angular dependence on the magnetization orientation of the ferromagnetic layer 
is substantiated by a phenomenological theory based on the symmetry of the effective spin-orbit 
interaction induced by the topology of the spiral magnetic order in the multiferroic barrier. 
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With the recent advances in fabrication, patterning and measurement techniques of arti- 
ficial nanostructures, thermoelectricity has gained renewed interest with a particular focus 
on potential applications [l 3 j , specially for the development of electronic devices with low 
power consumption. Thermoelectric transport is well studied theoretically and experimen- 
tally for a variety of systems such as quantum dots , molecular junctions |6||, and mag- 
netic tunnel junctions [v-9| . Recently, a qualitatively original phenomena, the Spin-Seebeck 



effect, was discovered by Uchida et al. 
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In a ferromagnetic material (NigiFeig) a 



heat current leads to a pure spin current without a charge current. This surprising effect 
has by now been confirmed experimentally for insulating ferrimagnets (LaY 2 FesOi2 in 12)). 
and for ferromagnetic semiconductors (GaMnAs [13|) as well. In these experiments, the 
spin-Seebeck coefficient was found to be quite smaller than the charge- Seebeck coefficient. 
However, in contrast to bulk ferromagnets, the spin thermopower in a quantum dot is theo- 
retically predicted to be on the same order of magnitude as the charge thermopower 
Along this line we inspect in this work the thermoelectricity in a normal-metal/helical- 
multiferroic/ferromagnetic heterojunctions such as Pt/TbMnOs/SrRuOa (cf. FigfTJ. The 
coexistence of coupled electric and magnetic order parameters in multiferroics [l^] holds the 



promise of futuristic opportunities for spintronics devices |15l . [16] with the particular advan- 
tage of being electrically 17) and/or magnetically 18) controllable. The essential point is 
that due to the topology of the local magnetic moments in multiferroic insulator, a travers- 
ing carrier experiences an effective spin-orbit interaction [16]. As shown by Hatimi et al. 
the thermoelectric effects significantly depend on the relative angle of the magnetization in 
neighboring magnetic layers of the ferromagnetic /normal- metal/ferromagnetic heterojunc- 
tions. Moreover, in the presence of the spin-orbit interaction, we expect that the tunnel 
magnetoresistance should show a spatial anisotropy in a ferromagnetic/insulator/normal- 
metal hetero junction 16, Q, 0]. Similarly, we find that the non-collinear magnetic order 
in the oxide together with the induced spin-orbit coupling result in uniaxial spin and charge 
thermopower with a C%v symmetry. 

A sketch of the system is shown in Fig{TJ A magnetic tunnel junction consisting of 
an ultrathin helical-multiferroic (MF) barrier (usually an oxide) is sandwiched between a 
normal metal (NM) lead and a ferromagnetic (FM) conductor. No voltage is applied across 
the junction. A steady temperature difference (T NM ^ T FM ) between the FM and NM 
electrodes is induced. As the transmission is spin dependent, the temperature gradient may 
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lead to a spin accumulation in the electrodes, which generally results in a nonzero spin- 
voltage bias AV S = (6}J,fm — 8fasrM)/e, where e is the electron charge and 5 fa = Sfa-f — 5fa± 
with 5 fan being the electrochemical potential of the spin a to the right or to the left of the 
oxide barrier. In the linear response regime (which is assumed here), the spin dependent 
current is written as [7J, 

I a = eL 0a Afa + ^L la AT (1) 

where Afar = eAV a is the difference in the chemical potentials of the two leads in the spin 
channel cr = ±l(or f4-)> an d AV a = AV e + aAV s with AV e being the charge bias. AT is the 



applied temperature difference. The kinetic coefficients L ncr are defined as 
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t 1 l dE (f urr(F) df{E) (?) 

~ ~h J 27r _ ~dE~ (2) 

where f{E) is the Fermi-Dirac distribution function. T a (E) is the spin-dependent trans- 
mission probability through the tunnel junction. The charge and the spin current are thus 
given as I e = 1^ + J^, I s — 1^ — 1^, respectively. We assume the spin-orbit interaction in 
the electrodes to be negligible, hence the two spin channels are independent. A nonzero 
spin current can be driven by the temperature difference through the system in absence of 
a charge current I e = as observed in the spin-Seebeck effect 10]- In the present study, 
we derive the thermoelectric and thermospin coefficients in the presence of the spin accu- 
mulation. To find the spin-dependent Seebeck coefficient 5*0-, we inspect the limit of the 
simultaneous vanishing of both the spin current and the charge current; or equivalently a 
zero charge current in each spin channel, i.e. I a = 0, 

_ AT _ eT Lqu ' [6) 
Analogously, the spin thermopower S s and the charge thermopower are calculated as, 



S, = ^ = l{S t + Si). (5) 

Neglecting the relatively much smaller potential modification due to the depolarizing field 
in the multiferroic barrier, and assuming that the barrier potential has a rectangular shape 
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with the height Vq, the Hamiltonian describing the tunneling across the hetero junction reads 



16|, H = H m + Hmf ■ H m stands for the itinerant carriers in the two metal electrodes, 



h 2 

H m = V 2 -Q(z- d) Am <x (6) 

2m e 

where cr is the vector of Pauli matrices, m = [cos 0, sin <ft, 0] is a unit vector defining the in- 
plane magnetization direction in the ferromagnet with respect to the [100] crystallographic 
direction, and A describes the Zeeman splitting in the FM electrode. Q(z) is the Heaviside 
step function. m e is the free-electron mass and d is the thickness of the barrier. In the oxide 
insulator, the carrier dynamics is governed by the exchange model, 

h 2 

H MF = V 2 + Jn r ■ cr + Vq, for < z < d (7) 

2m 

where m is the effective electron mass of the oxide ( m/m e ~ 3). Jn r is the exchange field, 



where n r is given by the multiferroic oxide 
by the integer number /) along the z-axis 



ocal magnetization at each spiral layer (labeled 
22], i.e., n r = (— l) z [sin 9 r , 0, cos 6V] with 9 r = q m r 
and q m = [q, 0, 0] being the spiral spin-wave vector. In effect the exchange coupling acts on 
the electron as a non-homogenous magnetic field. Performing a local unitary transformation 
within the barrier jisj, we conclude that the influence of the barrier amounts to the spin- 
dependent potential 

H e s f J = w(6 r ,k) ■ cr (8) 

where 

w(0 n k) = [J(z) sin 8 r , qk x , J(z) cos 6 r ] (9) 

and q = ^-Q- This effective spin-orbit interaction results in a tunneling anisotropic magne- 
toresistance (TAMR) effect [id]. Hence, we can expect a similar anisotropic behavior of the 
thermopower, as well. 

In the present study, we assume a barrier of two to five layers (as in j^]) such that the 
effective spin-orbit interaction HIjf throughout the multiferroic barrier is reduced to the 
plane of the barrier, H^f = w(0 r , k) • a5(z) with w(0 r ,k) = [J sm9 r ,qk x , Jcos9 r }. J and 
q are renormalized exchange and resonant spin-orbit coupling parameters, q ~ qdm e /m and 
J ~ (J(z))d referring to space and momentum averages with respect to the unperturbed 
states at the Fermi energy. In the following, we treat J and q as adjustable parameters. The 
transmissivity of a spin-cr electron through the multiferroic tunnel junctions reads 



T a (E,k h 9 r ) = M 



h h- 

— \+ I 2 _|_ —\i _| 

K K 



(10) 
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where the transmission (t a>cr and t a ^) coefficients can be analytically obtained by solving for 



the scattering states in the different regions 
in NM and FM subsystems, respectively, 



16]. k and k a are the transverse wave vectors 



K = d(E + aA)/£--kl (12) 

where ky denotes the conserved electron momentum parallel to the junction interfaces. 

Introducing the transmissivity T a (E, ky, 9) into Eq.([2]), the kinetic coefficients are rewrit- 
ten as, 

^<*> = -I I dE W?T. (z-^UEM.ef-m (13) 



Based on a general symmetry considerations of the spin-orbit interaction 24| and phenomeno 



logical calculations [16|] the angular-dependence of L ncr (0) is found to exhibit a two fold sym- 
metry, ~ cos 2(f). Consequently, the charge and spin thermopower are spatially anisotropic. 

We performed numerical calculations for k^T = 4meV with [i = 5.0eV, A = 2eV, 
Vo = 0.5eV, and d = lnm. Figj2] presents the dependence of the thermopower on the 
magnetization direction in the FM electrode. As evident from the numerical results, the spin 
and the charge thermopower show the symmetry, as follows from the phenomenological 
model. We note, the spin thermopower S s is about three orders of magnitude smaller than 
the charge thermopower ^e, which is different from the case of a quantum dot where S s 
can be as large as S e {4]. S s changes sign as voltage induced in the minority spin channel 
is higher. For the charge thermopower, the amplitude of the angular-dependence of S e 
is quite small, which is on the same order as the tunnel anisotropic magnetoresistance in 



Fe/GaAs/Au tunnel junction that have been recently realized experimentally [25|. However, 
S s changes with from positive to negative, we have a quite large tunnel anisotropic spin 
thermopower, [(S s ((j)) - S s (0))/(S s ((j)) + S s (0))} max ps 18. 

Summarizing, we studied the angular-dependence of the thermoelectric transport through 
the helical-multiferroic tunnel junctions, both the spin and charge thermopower are found to 
exhibit an anisotropic behavior due to the spiral magnetic order together with an induced 
spin-orbit interaction in the multiferroic spacer. Based on the magnetoelectric coupling, 
the strength of the effective spin-orbit coupling is electrically/magnetically controllable and 



thus the spin and charge thermopower in the helimagnetic tunnel junctions. For practical 
applications a multilayer configuration might be more appropriate to enhance the effect. 
This work is supported by the German Science Foundation, DFG through SFB762 -B7- 

functionality of oxide interfaces. 
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FIG. 1. (Color online) Schematic of a magnetic tunnel junction composed of a multiferroic oxide 
insulator barrier with spiral magnetic ordering (white arrows), a normal metallic (NM) lead and 
ferromagnetic (FM) electrode at different temperatures, Tnm 7^ Tfm- The vector m indicates the 
magnetization orientation specified by the angle <j) in xy (FM) plane with respected to the reference 
crystallographic axis (x-axis). The zx plane refers to the spiral plane of a multiferroic oxide. 




(/> 

FIG. 2. (Color online) Spin and charge thermopower as a function of the magnetization orientation 
cf) in FM layer at k^T = 4meV. Other parameters are chosen as /i = 5.0eV , A = 2eV, Vo = 0.5eV , 
d = lnm, J = leV and q = jj- with a = 5Abeing the lattice constant of the oxide. 
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